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Abstract. We consider the Klein- Gordon equation on two half-axes connected at their origins. 
We add a potential that is constant but different on each branch. In a previous paper, we studied 
the L°°-time decay via Hormander's version of the stationary phase method. Here we apply 
these results to show that for initial conditions in an energy band above the threshold of the 
tunnel effect a fixed portion of the energy propagates between group lines. Further we consider 
the situation that the potential difference tends to infinity while the energy band of the initial 
condition is shifted upwards such that the particle stays above the threshold of the tunnel effect. 
We show that the total transmitted energy as well as the portion between the group lines tend 

]/2 

to zero like a 2 in the branch with the higher potential ai as 02 tends to infinity. At the same 
time the cone formed by the group lines inclines to the i-axis while its aperture tends to zero. 



1. Introduction 

In this paper we study the energy flow of waves in two coupled one-dimensional semi-infinite 
media having different dispersion properties. Results in experimental physics 0[9], theoretical 
physics [7] and functional analysis [H E] describe phenomena created in this situation by the 
dynamics of the tunnel effect: the delayed reflection and advanced transmission near nodes 
issuing two branches. Our purpose is to describe the influence of the height of a potential step 
on the energy flow of wave packets above the threshold of tunnel effect. 

We consider the following setting: let N±, N2 be disjoint copies of (0, +00). Consider numbers 
0^ satisfying < a\ < ci2 < +00. Find a vector (ui,U2) of functions u k ■ [0, +00) x N k — > C 
satisfying the Klein-Gordon equations 

[8 2 -d 2 x + a k ]u k {t,x) = 0, k = 1,2, 

on N\ , N2 coupled at zero by usual Kirchhoff conditions and complemented with initial conditions 
for the functions u k and their derivatives. 

Reformulating this as an abstract Cauchy problem, one is confronted with the self-adjoint 
operator A = (—d 2 + a\, —d 2 + 02) in L 2 (Ni) x L 2 (N2), with a domain that incorporates the 
Kirchhoff transmission conditions at zero. For an exact definition of A, we refer to Section [2j 
The problem described above can be reformulated as 

u(t) + Au(t) = 0, 

u{t) e D(A), (l) 

for all t > together with initial conditions. It is well known that the following expression is 
invariant with respect to time for solutions of ([T|): 



E(u(t,-)) = ^\\u(t,-)\\ 2 H + (Au,u) H y 
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In Section [2] we recall the solution formula that was proved in [2] by an expansion in generalized 
eigenfunctions in the more general setting of a star shaped network with semi-infinite branches. 

In Section [3] we recall our result on L°°-time decay proved in [3j. There we obtained the 
exact L°°-time decay rate c • t~ l l 2 in the group velocity cones together with an expression for 
the coefficient c for initial conditions in a compact energy band in (ci2, 00). In the present work 
we refine an estimation from below in this context. 

In Section 2] we use the preceding results to estimate the L 2 -norm of the outgoing solution on 
N2, which is a part of the total energy given in ([2]). We suppose that the initial condition belongs 
to an energy band above the threshold of the tunnel effect, so that the solution propagates in 
the branch with the higher potential A 2 . We consider the L 2 -norm of the solution at time t 
both on the whole branch N2 and inside the cone delimited by the group lines corresponding to 
the bounds of the energy band. It turns out that the first norm has an upper asymptotic bound 

— 1/2 

and and the second one an upper and lower asymptotic bound which behave as o 2 , if the 
height of the potential step 02 tends to infinity. This implies that the ratio of the energy on the 
whole branch and the energy between the group lines is time asymptotically confined in a finite 
interval above 1 which is independent and a%. 

These results might be interpreted in terms of quantum mechanics as follows: a relativistic, 
massive particle without spin in a one-dimensional world is submitted to a potential step at 
the origin. It is supposed to have enough kinetic energy to overcome the step with a fixed 
remaining energy. Classically, the particle should leave the potential step with a velocity which 
is independent of the height of the step. Our results show that in the quantum mechanical model, 
when the height of the potential step tends to infinity, the velocity of the outgoing component 
of the particle tends to zero while the particle is more and more localized. At the same time 
the total outgoing energy tends to zero while its ratio to the energy inside the cone delimited 
by the group lines remains time asymptotically in a constant finite interval above one. Similar 
estimates should be possible for the other parts of the total energy ([2]). 

The last observation might lead to the idea of viewing the total energy as being subdivided in 
a part inside the group line cones and a part outside these cones. The ratio of these two energies 
is quite independent of the experimental configuration, but depends mainly on the chosen energy 
band. 

The results of this paper are related to results in experimental physics, theoretical physics 
and functional analysis (spectral theory, asymptotic estimates, analysis on networks, cf. [3] and 
the references cited there). For example in [2] we obtained information on the splitting of the 
energy flow near zero. In [1] a (not optimal) estimation for the L°°-time decay rate has been 
obtained but without any information on the localization of the energy. In [5] the relation 
of the eigenvalues of the Laplacian in an L°°-setting on infinite, locally finite networks to the 
adjacency operator of the network is studied. In [10], the authors consider general networks 
with semi-infinite ends. They give a construction to compute some generalized eigenfunctions 
but no attempt is made to construct explicit inversion formulas. 

2. A SOLUTION FORMULA 

The aim of this section is to recall the tools we used in [2] as well as the solution formula of the 
same paper for a special initial condition and to adapt this formula for the use of the stationary 
phase method in the next section. 

Definition 2.1 (Functional analytic framework). 

i) Let Ai, A2 be identified with (0, +00). Put N := Ni x A 2 , identifying the endpoints 0. 



ENERGY FLOW ABOVE THE THRESHOLD OF TUNNEL EFFECT 3 

ii) Two transmission conditions are introduced: 

(T ): (u!,u 2 ) G C(7h) x C(Jh) satisfies ui(0) = u 2 (0). 

This condition in particular implies that (ui,u 2 ) may be viewed as a well-defined function 
on N. 

2 

(Ti): (ui,u 2 ) G C7 1 (iVi) x C 1 ^) satisfies ^9^(0+) = 0. 

k=l 

iii) Define the Hilbert space H = L 2 (N\) x L 2 (N 2 ) with scalar product 

2 

(U,V) H = ^2(uk,V k )L2(N k ) 
k=l 

and the operator A : D(A) — > H by 

D{A) = |(ni,n 2 ) G H 2 (N{) x H 2 {N 2 ) : («i,« 2 ) satisfies (T ) and (Ti)}, 
^(ni, u 2 ) = ^2^2) = (y-dlu k + 1 



Ofc-Ufc 



ft=l,2 



Note that, if a% = a 2 = 0, then A is the Laplacian in the sense of the existing literature, cf. 



Definition 2.2 (Fourier-type transform V). 
i) For k =1,2 and A G C let 



6(A) £i(A) 



£ fe (A) := \/\-a k as well as sj(A) := -7-7TT and s 2 (A) :- . 

SI (A) £2 (A) 

Here, and in all what follows, the complex square root is chosen in such a way that 

Vr • = ^Jre % ^l 2 with r > and G [— 7r, 7t). 

ii) For A G C and j,k G {1,2}, we define generalized eigenfunctions i 7 ^'*' : N — > C of A by 
F A ±J (x) :=F±' k j (x) with 

f F\k{x) = cos(£j(A)x) ± isj(A) sin (£,,■( A) x), for A; = j, 

I F x,k( x ) = exp(±i£ fc (A)x), for /c / j. 

for x 6 iVfc. 

iii) For I = 1,2 let 

(O, ifA<^, 
«W := S ft(A) if n. \ 

iv) Considering and (72 as weights for our L 2 -spaces, we set L 2 := L 2 ((ai, +00), q\) x 
T 2 ((a2, +00), q 2 ). The corresponding scalar product is 



2 /• 

(F, G) q := £ / %(A)F fc (A)G fc (A) d\ 

i. 1 J (at ,+00 1 



k=1 J(a k ,+oo) 

:=(F,F)l /2 . 

v) For all / G L l (N, C) we define Vf : [a x , +00) x [a 2 , +00) -> C by 



1/2 

and its associated norm \F\ q := (F,F) q / 



(Vfh(X):= / /(x)(F-' fc )(x) dx, fc = l,2. 
./TV 
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In [2], we show that V diagonalizes A and we determine a metric setting in which it is 
an isometry. Let us recall these useful properties of V as well as the fact that the property 
u £ D{A ) ) can be characterized in terms of the decay rate of the components of Vu. 

Theorem 2.3. Endow C c °°(iVi) x C™(N 2 ) with the norm of H = L 2 {Ni) x L 2 (N 2 ). Then 

i) V : C^°(Ni) x C^°(N 2 ) — > L 2 is isometric and can be extended to an isometry V : H — > 
L 2 , which we shall again denote by V in the following. 

ii) V : H — > L 2 is a spectral representation of H with respect to A. In particular, V is 
surjective. 

iii) The spectrum of the operator A is o~(A) = [a%, +oo). 

iv) For I £ N the following statements are equivalent: 

(a) u G D(A l ), 

(b) A i->- X l (Vu)(X) £ L 2 , 

(c) X l (Vu) k (X) £ L 2 ((a k ,+oo),q k ), k = 1,2. 
We are now interested in the Abstract Cauchy Problem 

(ACP) : u a (t) + Au(t) = 0, t > 0, with u(0) = u , u t (0) = 0. 

Here, the zero initial condition for the velocity is just chosen for simplicity, as we will not deal 
with the general case in this contribution. 

By the surjectivity of V (cf. Theorem l2.3l (ii)) there exists an initial condition uq G H satisfying 

Condition (A): (Vu ) 2 = and (Vu )i G C 2 ((a 2 , oo)). 

Remark 2.4. i) For -uo satisfying (^4) there exist a 2 < A m j n < A max < oo such that 

supp(Vuo)i C [A min ,A max ]. 

ii) If u G H satisfies (A), then u G D{A°°) = f] D(A l ), due to Theorem O (iv) , since 

i>i 

X i y X l (Vu) m (X) G L 2 ((a m , +oo), q m ), m = 1,2 for all I G N by the compactness of 
supp(V« ) 

Theorem 2.5 (Solution formula of (ACP) in a special case). Suppose that uq satisfies Condition 
(A). Then there exists a unique solution u of (ACP) with u G C l ([0, +oo), D(A m / 2 )) for all 
I, m G N. For x £ N 2 we have the representation 

u 2 {t,x) = -(u + (t,x) +u-(t,x)) 

with 

u ± (t,x):= / e ±i ^ i (7l (A)e-^^(yn ) 1 (A)dA. (3) 

J ^min 

Proof. Since vq = ut(0) = 0, we have for the solution of (ACP) the representation 

u(t) = V~ l cos(\/~Xt)Vu Q . 

(cf. for example [U Theorem 5.1]). The expression for V' 1 given in [2J yields the formula for 
u±. ' □ 

Remark 2.6. A solution formula for arbitrary initial conditions which is valid on all branches is 
available in [2]. This general expression is not needed in the following. 
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3. L°°-TIME DECAY 

Next, we quote the result on L°°-time decay of solutions to the problem (ACP) from [3]. Here, 
we only consider special initial conditions that are localized in energy in a compact interval 
contained in (02,02 + !)• For these it is possible to give very explicit estimates for all the 
constants that appear in a asymptotic expansion of the solution to the order t 1 / 2 . 

Theorem 3.1. Let < a < j3 < 1 and ip £ C 2 ((a,/3)) with \\ip\\oo = 1 be given. Setting 
:= ip(\ — a 2 ), we choose the initial condition uq £ H satisfying (1^0)2 = and (Vuq)i = t/j. 
Furthermore, let ii + be defined as in Theorem \2.b\ 

Then there is a constant C(vp,a,f3) independent of a\ and a 2 , such that for all t £ M + and 
all x £ N 2 with 



we have 

\u + (t,x) - H(t,x,u ) ■ t~ 1/2 \ < C(ip,a,P) ■ r 4 , 

where 

H(t,x,u ) :=e- i ^ > t > x \2m)y 2 4 /4 h^x) h 2 (t,x) {Vu )i(a 2 + p 2 ) 

with 



ip(p, t, x) := s/a 2 +p 2 t - px, po 



a 2 x 2 
(t 2 - x 2 ) ' 



. ( . , ( (t/x) 2 . , V(a 2 -ai)((t/x) 2 -l) + oi 
hi(t,x) := I 7~~2 7 , h 2 (t,x) := — -3. 

\ {t/x > -V (V(«2-ai)((t/x) 2 -l) + a 2 + ^) 

It holds further 

TTU \\ * I R\ r^T V / ?( Q 2+/3) 3/4 /r— r -1/4 

H [t, x, uo)\ < g{a\,a 2 , p) := y Itt ~ \>2Trp a 2 as a 2 —> +00. 



Proof. This is contained in Equation (8) from the proof of Theorem 3.2 and in Theorem 4.1 in 

0- □ 

Remark 3.2. i) Note that (HJ) is equivalent to 

1>min < V(t,x) := (t/x) 2 - 1 < U ma x • 

, a 2 a 2 a 2 a 2 

where u min := = — and v max := = — . 

Amffl — 0, 2 p A m in — 02 « 

ii) For later use we define p min := £ 2 (A min ) and p max := &(\ ma x)- 



4. Energy flow 

In this section we use the asymptotic expansion from section to estimate the outgoing 
solution from above and below in the cones given by the group velocities corresponding to the 
bounds of the energy band of the initial condition. This leads to time independent asymptotic 
upper (Theorem 14. 4p and lower (Theorem 14. ip estimates of the L 2 -norm of the solution on the 
space interval inside the cones. 

Further an upper estimate of the L 2 -norm of the outgoing solution on the whole branch N 2 
is obtained using PlanchereFs theorem (Theorem 14. 2p . 
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These informations lead to our main result (Theorem I4.3p . where we give an upper bound 
for the ratio of the L 2 -norm on the whole branch N 2 and the L 2 -norm inside the cone, which is 
asymptotically independent of the height of the potential step. 

Theorem 4.1. In the setting of Theorem \3.1\ suppose that > m > for fj, £ [a' , j3'] with 
a < a' < (3' < (3. Then we have 

i) the lower estimate for the coefficient oft~ x l 2 : 



\H(t,x,u )\ > f(a 2 ,a,(3) ■ m 

:=V^af (£ + l) 3/4 1 /02 " ai 

~ V2ira a 2 m, as a 2 
for all (t, x) satisfying 



\fa~2 V a 



+ 1 



-2 



a 2 — ai 
a 



+ 1 + 1 m. 



00 



a 2 + P' t ^ ja 2 + a' 



y P' x V a' 
ii) and the lower estimate for the solution 

Ve > 3t > Vi > t \u+(t,x)\ > (/(a 2 , a, /3)m - e) t~ 1/2 
for all (t, x) satisfying 



a 2 + /?' t ja 2 + a 1 



ft' 



a' 



Proof. Note that it is always possible to choose the initial condition in the indicated way, thanks 

to the surjectivity of V, cf. Theorem 12.31 ii). 

(i): 

Theorem 13.11 implies 



\H(t,x,uo)\ = V2Tra2^hi(t, x)\h 2 (t,x)\ ■ [|(Vuo)i| 



We estimate 



\h 2 (t,x) 



y/{a 2 - ai )((t/x) 2 -l) + a 2 



> 



(\A°2 - ai){(t/x) 2 - 1) + a 2 + 
V( a 2 ~ ai)v max + a 2 



> 



(\/(a2 - OlHnax + 02 + v^) ' 
V(«2-ai)^ + a 2 



(V(«2- ai)^ + a 2 + ^i) 



+ 1 



\Z&2 ( / Q2-Q1 



+ 1 + 1 



~ as a 2 — )• 00 

a2 
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Here we used, that b i-> nq^p is decreasing for 6 > c > 0. Further 

( Amax \ 3/4 n 

, Amax ~ a2 J = (— + 1) 3/4 -> 1 as a 2 -»• oo 
7 Amax 1 / a 2 

This implies (i). 
(ii): 

Using the lower triangular inequality we find that Ve > 3to > Vi > to : 

> |it+(i,x) - H(t,x,u ) ■ t~ 1/2 + H(t,x,u ) ■ r 1/2 \ 

> \H(t,x,u )\ ■ t~ l/2 - \u + (t,x) - H(t,x,u ) ■ r l/2 \ 

> (f(a 2 ,a,P)m-s) r l/2 

for all (t, x) satisfying 



'02 + < t < /a 2 + «' 



This shows (ii) □ 
Theorem 4.2. Suppose the setting of Theorem VJ.1\ Then we have 

\\u+(t, -)\\v>r Na ) < = =-r=||^|| I ,2 ( ( Q ,^)) ~ ^2||V'llL2(( a ,/3))a^ 1/2 ^ a 2 -»• oo 

Proof. In the expression ([3]) for u + we substitute p := £ 2 (A). This yields 

/Pmax > 
e^ S2+p2t qi {a 2 + p 2 )e-^(yu )i(a 2 + p 2 )p dp. 
i min 

Interpreting this integral as a Fourier transform we find by the Plancherel theorem 

II«+(*,0IIl2(jv 2 ) < IK(V)IIl2(R) = ||p^poi(a 2 + p 2 )(Vu )i(a 2 +P 2 )X[p min , Pmax ](p)|| L 2 (K) 

Now, we use that supp((Uuo)i) is contained in the interval [a 2 + a,a 2 + 0\, so only the range 
\[ol < P < v 7 /^ is relevant. For these values of p we find 



pqi(a 2 +p ) =p — < - = < 



(a/o^ - ai +p 2 +p) 2 x/az - ai +p 2 V«2 - ai + a 
So, we have 

\\u+(t,-)\\ L 2 {N2) < ||(Vno)i(o2 +P 2 )\\ L2((v . _ » 

and substituting back A = a 2 + p 2 we obtain 

11(^0)1(^+^)11^ = / |(yu ) 1 (A)| 2 ^_ = / |,/,(A)| 2 ^_. 

^ UPmm,PmaxjJ A min y/X-a 2 Ja 2 +a V A - a 2 

Setting, finally, /i = A — a 2 we end up with 

ik(v)iu w < r J2== ( f\m?^=f' 2 < , = u\\l^)) □ 

V fl 2 - oi + ol \J a ■y/JiJ yja 2 — ai + a\?a u 

Theorem 4.3. In i/ie setting of Theorem \3.1\ suppose that ip(p) > m > /or /i £ [c/,/3'] 
a < a' < P' < (3. Then we have 
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i) the estimate from below: Ve > 3to > Vt > to 



I M*. •) > (/K a, - E ) 2 (^^7 " + a , 



(Y 



where 



I'f-- 



o 



t. 



«2 + a' V °2 + P' 



and f is defined in Theorem \3.1\ (i). 
ii) the estimate for the ratio of the global and local energy: 
Ve > 3t > Vt > t 



\\u + (t, 



\\u + (t, ■ 
iii) and 



\l 2 (n 2 ) . \fP \\il>\\i?{{fx,p ))(., m N_1 
I < / = T7=[f(a2,a,p)m-e 



/3' 



a2 + /3' V a2 + a' 



a' \-l/2 



lim sup 

t— >oo 



ll«+(t, 



< 



1 oil^j/) 



a/oi — a2 + a^/a 



f(a 2 ,a,ft)m 



ft' 



0.2 + ft' V a 2 + a' 



a' \ -1/2 



(2^)- 1 /2 m -i /3 i/2 a -3/4 (v ^7_ ^)-^ 2 U\\ LH{aM as 



Q2 — > OO 



Proof, (i): 
Follows from 



||«+(t, 



lL2(/ t ') 



n > ( inf \u+(t, ■)\) • |l t '| 



and Theorem 14.11 (ii). 

(ii) : 

Follows from (i) and Theorem 14.21 

(iii) : 

Direct consequence of (ii) and Theorem 14.11 (i). 

Theorem 4.4. Consider the setting of Theorem \3.1\ 

i) Ve > 3t x > Vt > h : 

\u + (t,x)\ < (g(a 1 ,a 2 ,/3) + e) r 1 ' 2 

for all (t, x) satisfying 



1 0,2+ ft t ja2 + a 



ft ~ X 



a 



ii) Ve > 3ti > Vt > t x : 



I M*. •) < (9(01,02, ft) + ef -\la 2 + a 



a 



□ 
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limlof||u + (t,-)||^ f j» < g(ai,a 2 ,f3) 2 (J — ^-5- J 

~ 2-7r/3(y / /3 — yfa) 1 as 02 00 

Proof, (i): 

TheorenTlO implies that Ve > 3t x > Vt > *i : 

\u + (t,x)\ < \u+(t,x) - H(t,x,u )r 1/2 + H(t,x,uo)t~ 1/2 \ 

< C(^,a,/3)t" 1 + \H(t,x,u )\t- 1/2 

< (g( ai ,a 2 ,P) +e) r 1 / 2 

for (t, x) in the cone indicated there. 

(ii) : 

Follows from estimating the square of the L 2 -norm against the square of the maximum of the 
function (using (i)) times the length of the integration interval. 

(iii) : 

Direct consequence of (ii). 

□ 
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